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^ Abstract 

■ We extend the predictive soft-wall AdS/QCD model to a thermodynamic model by considering 
the black hole metric. A modified bulk vacuum expectation value and a modified coefficient for 
the quartic term in the bulk scalar potential are introduced to obtain a smooth dilaton solution for 
computing the spectral functions of the vector and axial- vector mesons. It is demonstrated that the 

■ peaks appearing in the spectral functions characterize the thermal mass spectra of vector and axial- 
vector mesons, where the location of the peak moves to a lower value and the width of the peak 
becomes wider when increasing the temperature. We observe that the peak disappears completely 

r"! ! at the critical temperature around T^. = 200 MeV, which implies the deconfinement of quark and 

the restoration of chiral symmetry breaking. A numerical study by fitting the spectral function 
in terms of the Breit-Wigner form has been made to show how the peak dissolves quantitatively 
when the temperature is increased to the critical point. 
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I. INTRODUCTION 



The reproduction of the extremely early moment after the big-bang is one of the biggest 
purposes and achievements in the present high energy physics. The experimental exploration 
is expected to be carried out by the groups at LHC and RHIC, as well as many other groups 
including FAIR, J-PARC, JLab and RIBF. It will involve with the following basic problems: 
the novel matter phase of quark-gluon plasma, confinement of colors and quarks, chiral 
symmetry breaking and restoration, unstable nuclides, inner structure of hadrons and the 
force between hadrons, color glass condensation and superconductor. Unlike the perturbative 
QCD which is well understood based on the property of asymptotic freedom of QCD[i|, 0], 
all the studies will lead us to understand deeply the properties of nonperturbative QCD 
as they all are relevant to the strongly coupled QCD. So far, the strongly coupled QCD 
is mainly described by the chiral dynamical models of QCD with dynamically spontaneous 
symmetry breaking jsHsl, lattice QCD and holographic QCD (or AdS/QCD)j6[. 

The holographic QCD has been developed based on the AdS/CFT duality to character- 
ize the strongly coupled QCD. AdS/CFT was initiated as the duality between the weakly 
coupled type IIB supergravity on AdS^ x space-time and the strongly coupled D = 4 
A/" = 4 super Yang-Mills theory 0]. It is known that D = 4 Af = 4 super Yang-Mills the- 
ory has conformal symmetry and no field associated with fundamental representation. Thus 
AdS/CFT itself cannot correspond to the realistic QCD. The attempt for making AdS/CFT 
approach to the realistic QCD was first pursued in ref.j6|, where the study has been paid to 
construct phenomenologically the model in searching for the bulk gravity matching with the 
realistic QCD. Such kind of bottom-up approach is usually called as holographic QCD or 
AdS/QCD models which include hard-wall AdS/QCD model and soft- wall AdS/QCD 

model [l2-T^. It is also interesting to observe the correspondence between matrix elements 



obtained in AdS/CFT with the corresponding formula using the light-front representation as 



shown in Refs. |17l4l9l|. It may distinguish from the top-down model resulted from solutions 
of string theories like D3/D7 |20j, D4/D8/D8 |2l| and D4/D6/D6 Q. 

In AdS / QCD models, the cutoff is put artificially around the horizon in the radial coor- 
dinate of the bulk gravity. The location of such a cutoff corresponds to the inverse of QCD 
scale 1/Aqcd in the dual field theories. Here the radial coordinate in the bulk corresponds 
to the energy scale in the dual field theories, thus the cutoff has a role as the IR cutoff in 
the dual field theories. There is a difference in how the cutoff is put between hard-wall and 
soft-wall AdS/QCD model. In the hard-wall model, the range of the radial coordinate is 
restricted sharply at the cutoff in the integration of Lagrangian. Whereas in the soft-wall 
model, the action begin to be suppressed gradually from around 1/Aqcd to the horizon by 
introducing the dilaton which is regarded as a background field. 

One of the main topics in the holographic AdS/QCD is to consistently incorporate both 
the dynamics of chiral symmetry and the linear dependence of squared mass of mesons on 
spin and excitation number, namely the dynamically spontaneous chiral symmetry breaking 
and the linear confinement of QCD. The hard-wall model can succeed in the realization of 
the linear dependence on lower excited states and lower spins, while the linear dependence 
deviates at higher excited states or higher spins j23[. The soft- wall model was motivated 



to improve this situation [12|, so that the linear dependence was found to be satisfied for 
higher excited states and spins, but the resulting chiral condensation is proportional to 
quark mass, which is inconsistent with the realistic QCD. To further improve the situation, 
a higher order quartic term was introduced in the potential of the bulk scalar field, which 
modified the proportional between chiral condensation and quark mass, but the bulk scalar 



potential field is not bounded [14 



It has been shown in ref . [15|, Il6| that by simply modifying the bulk gravity at the infrared 
(IR) region in the soft- wall AdS / QCD and keeping the conformal invariance to be unchanged 
at the UV region as required from the property of QCD, it enables us to build a predictive 
AdS/QCD model which can provide a consistent prediction for the mass spectra of all light 
mesons. As a consequence, the mass spectra for both the groundstate mesons and resonance 
mesons can match well with the experimental data. 

Another interesting topics in the holographic AdS/QCD are to investigate the finite 
temperature effects. In ref.jlS], in both of the hard- and soft-wall models, confine- 
ment /deconfinement transition temperatures were computed based on the evaluation of the 
free energies. Ref. j25| is a study of meson in the hard-wall model at finite temperature. 
The finite temperature effect on the spectrum of glueballs or mesons in the soft-wall model 



was studied in ref. |26l-l28l| . In Refj29[, a holographic model for charmonium was built in 
the soft- wall model. In refs. [isl, Soj , both mesons and gluons were investigated in the soft- 
wall model. In ref. , the spectral function was computed by taking the D3/D7 setup 
which is dual to TV = 4 Super Yang-Mills theory with fundamental matters at finite baryon 



density. In ref.jSJ], the spectral function was computed holographically by modifying the 



D3/D7 setup with finite baryon density |3ll-l33| such that the dual field theory can have 
electric background. The setup with D4 and D4-branes on AdSg was considered in ref. [35| . 
This is a non-critical string setup, which is considered not to have the problems originated 
from in critical string setups like DA/DA/DS. In such a circumstance, the spectral function 
was computed (hydrodynamics was also performed). 



In this paper, we are going to extend the predictive soft- wall AdS/QCD model [15|, [16 
to a thermodynamic model at finite temperature. In particular, we will examine finite 
temperature effect on mass spectra of vector and axial-vector mesons from the computation 
of the spectral function. The paper is organized as follows: In Secjlll we include the finite 
temperature effect into the predictive soft- wall AdS/QCD model in the usual way by simply 
considering a black hole metric. In sec JIIIt we present a numerical computation for the 
spectral functions of the vector and axial-vector mesons at finite temperature. It is seen 
that the peaks of the spectral functions at low temperature will appear like a spike, the 
locations of the peak correspond to the groundstate mesons and resonances, which are 
found to match with the ones at zero temperature. We will show how the peaks dissolve as 
temperature increasing, which then indicates the chiral symmetry restoration. In Sec JIVj we 
will qualitatively investigate the dissolution of the spectral functions by fitting to numerical 
results in terms of the Breit-Wigner form. The finite momentum effect in the spectral 
function is discussed in SecJVl At the end of the section, we make remarks on the options 
in parameters and profiles for the high order scalar interactions as well as the boundary 
behavior of the scalar VEV. Our conclusions and remarks are given in Sec JVIi Some formula 
used in our computation of spectral function are derived and presented inlAl 



II. PREDICTIVE ADS/QCD MODEL AT FINITE TEMPERATURE 



In this section, we will describe the extension of the predictive AdS/QCD model 
a thermodynamic model by simply considering the black hole metric [261 



11,0 



to 



It then turns 



out that if taking the previous bulk vacuum expectation value (bVEV) given in [15|, |16 



the resulting dilaton background will diverge on the horizon at finite temperature due to 
the modified bulk geometry in the IR-region and the quartic term in the scalar potential, 
thus a modified bVEV has to be introduced to avoid such a divergence in order to obtain 
a well-defined spectral function at finite temperature in the predictive soft- wall AdS/QCD 



model. 

Let us begin with considering the following black hole space-time as a bulk gravity: 



ds' = a\z)lf{z)dt'-Y.d^"-ir-^]^ (1) 



with the IR improved metric and the black hole metric 

a\z) = l/z' + fil and /(z) = 1 - (z/zo)'- (2) 

Here z has the relation with a black hole radial coordinate r via z = 1/r, it takes from 
zq = I/tq (inverse horizon) to (boundary). (In Sec JIIIl we will see that it is more convenient 
to use the variable u = z/zq.) As the boundary theories are independent of the AdS radius, 
we will set it as unity through this paper. The Hawking temperature of the black hole is 
given as T = l/{zo7r), which corresponds to the temperature in boundary theories. The 
mass scalae Uo characterizes QCD confinement and is fixed by the low energy physics of 



hadrons 15, 1 



We take the black hole space-time as the bulk geometry, despite what the boundary theory 
is in confinement phase, as in the field theory the mesons are considered as the bound states 
of strong interaction. In our treatment, the bulk geometry of the black hole space-time 
given in Eq. ([T]) is taken as the background. The dilaton is also taken as the background 
field, but it is obtained in the bulk geometry of the black hole space-time from solving the 
equation of motion for the scalar field with the required boundary conditions, the details 
can be seen below. Thus, our system can be a solution of equation of motion in the end. 
(While it is not known whether it is stable or metastable. One way for the confirmation of 



that is via the evaluation of the specific heat like ref.[28|. ) This is the different point from 
other studies 2^2^, where both of the dilaton and the black hole space-time are taken as 
the backgrounds for confinement phase. 

On the background geometry Eq.([T]), the considered field ingredients include the dilaton 
background the bulk scalar field X and gauge fields Al and of SUl{2) x SUr{2) 
gauge symmetry. Here X = XH" and Al^r = Alj^t" with T (a = 1,2,3) the SU{2) Lie 
algebra are as the dual operators in the boundary theories. The model is given as follows 
on the background Eq.([T]), 



5* — + 5'gauge) (3) 



where 

S. 



with 



- y"rf^Xv/^e-*(^)Tr(|A,Xp-\/(X)), (4) 

DmX = OmX + I {Al,mX - XAr^m) , (6) 
V{X) = m'M' + ^m'. (7) 



with rrix = —3. = 12tt L/Nc with Ai'c = 3 is a gauge coupling given in [12[. We use capital 



Latin and Greek character for the five dimensional and the four dimensional coordinates. 



respectively, M = xo,xi,X2,X3, z and /i = xo,xi,X2,X3. V{X) is the potential of the bulk 
scalar field. A is a coefficient for the quartic term. Here we will consider the case with or 
without the quartic term, i.e., A = and A 7^ 0. For A 7^ 0, its value at zero temperature 



was fitted from the meson mass spectra and found to be A = 9 [15|, Il6[ . 

Let us now discuss the bulk vacuum expectation value (bVEV) of the scalar field and the 
dilaton background at finite temperature. Taking the bVEV to be a function of the fifth 
dimension v{z) 

Xiz) = ^v{z) I2, (8) 

where I2 denotes 2x2 unit matrix. v{z) at zero temperature is given by Table HI HTl and UTTl 
Here we refer it as zero temperature part and denote it as Vq{z) below. 
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Ila 
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{B^Kl^dl'))'" 



TABLE L These are taken from 15|. Numerical values are given in Table HIl and lllli 



Model 


ruq (MeV) 


(j3 (MeV) (MeV) 


7 


la 


4.16 


275 


(2/^l)/3 


0.178 


lb 


4.64 


265 




0.136 


Ila 


4.44 


265 


(2/^d)/3 


0.153 


lib 


4.07 


272 


/^d/3 


0.112 



TABLE II: The numerical values in the case without the quartic term (Aq 
and = 445 MeV for all the models. 



0) m- C = \/3/(27r) 



Model mq (MeV) a 3 (MeV) (MeV) (MeV^ 



7 



Ila 
lib 



6.95 
6.79 



228 
229 



412 

548 



(2/i^)/3 
/"g/3 



0.30 
0.20 



TABLE III: The numerical values in the case with the quartic term (Aq 
for all the models. 



9) M. C = V3/(27r) 



At finite temperature characterized by the black hole metric, the dilaton field ^{z) has to 
be consistently determined from solving the equation of motion for any known bVEV v{z) 

a[z) f[z)v'[z) f{z)v'{z) V 27 



where the prime (') means a 2;- derivative. It is seen that there are terms with f{z) in the 
denominator in eq.(|9]), which leads to a divergent solution for the dilaton field ^{z) on the 
horizon {z — )■ 1) at finite temperature if taking the bVEV given in [l5| for zero temperature 
case. While the divergent solution of ^{z) will bring the difficulty in taking the in-falling 
boundary condition for the spectral functions. Note that the divergence appears only on 
the horizon at finite temperature due to the black hole metric f{z). At zero temperature 
= 1- 

To avoid such a divergence in the solution for the dilaton field ^{z), we may modify the 
bVEV v{z) to be 

v{z) = viiz)\nf{z) + voiz), (10) 



where vo{z) is taken to be the part at zero temperature and given in [15|. The part 



Vi{z) In f{z) vanishes at zero temperature and is introduced to obtain a finite solution for 
the dilaton field ^{z) by a proper vi{z). To fix Vi{z), we substitute the above modified 
bVEV into eq.dH]), which gives 

$'(^) = J'^Iil, l-a\z)v,iz) (ml + \\z)\ + f\z)v[{z)\ + ■ ■ ■ , 



vi{z)f'{z) 



here we write only the divergent term proportional to hif{z), and "■ ■ ■" denotes the finite 
terms at the horizon limit: z — )■ 1. Where we have defined the function J^{z) as 

:F{z) = -mla\z)v^{z) + nz)v[{z) (12) 

It is seen that for the case A = the solution for ^{z) becomes finite once J^{z) = 0. It is 
not difficulty to find from the condition J^{z) = that ^1(2;) should take the following form 



vi{z) = Cy exp 



(13) 



with c„ an integral constant. The above result holds for all the models discussed in 15|. We 
will make a detailed discussion below for the fixing of the constant c^. 

For the case A 7^ 0, there is an additional divergent term in Eg. (11 II) . When taking Eg. ( 1131) . 
the divergent term is given as 

$'(.) = -^^^a\z)v\z) + . . . . (14) 

its divergent order is seen to be (ln/(2;))'^. As a simple prescription, we may take A as the 
following form 

A = with p > 3, (15) 

l + CA(ln/(z))P ' ^ ^ 

with c\ a constant. Thus the resulting solution for $(2) becomes finite as long as p > 3. 
For simplicity, we take p = 3 in the present consideration. 

With the modified bVEV given in Eg. lfTOl) and the modified coefficient of the guartic 
term given in Eg. f lT5|l . we then arrive at a nonsingular solution for $'(2;) as shown in FigJTJ 




FIG. 1: The x- and y-axis mean u = z/zq and ^'[u). #'(u) is shown to be well-defined in the 
bulk entirely with the bVEV in Eq. (jlOp and the quartic coupling in Eq. (jl5p . T = 0.05 and 0.15 
GeV correspond to low and high temperatures in this paper. The correspondence between lines 
and parameters is given in Table ITVl Where the small windows give the magnified vicinity of the 
horizon, the red line (or thin solid line) and the blue line (or dashed line) show the divergence with 
behaviors as l/u^ and cubic logarithmic as expected from Eqs.Q and for = and c\ = 0. 
When Cv (green or dotted line) and c\ (magenta or thick solid line) are turned on, the divergences 
disappear. In the small windows, green line (dotted line) and magenta line (thick solid line) are 
overlapping each other. 



line 




c\ 


Ao 


red (thin solid) 











green (dotted) 


-0.1 








blue (dashed) 


-0.1 





100 


magenta (thick solid) 


-0.1 


-0.1 


100 



TABLE IV: The correspondence between the lines and the parameters used in Fig{TJ Other 
parameters are given in Table Hill The red line (thin solid line) and the blue line (dashed line) 
diverge as expected from Eqs.® and (fT^ . Where the green line (dotted line) and the magenta 
line (thick solid line) correspond to the modified bVEV given in Eq. (|10p and the modified coefficient 
given in Eq. (jl5p . 



Now let us further discuss the parameters c„ and c\. As they cannot be determined the- 
oretically in the present considerations, we are going to fix them from the phenomenological 
studies. Their numerical values will be given in the next section, here we shall first discuss 
the general considerations for how to fix them. 

It is seen that the parameter is introduced as the coefficient of the term vi{z) given in 



Eq.( !T3|) which is associated with the part of the VEV vi{z) In f{z) in Eq. (fTO!) . This part of 
the VEV plays the role of a regulator, otherwise the solution of dilaton will be divergent. 
Similarly, the parameter c\ is introduced in Eq. (fT5l) to modify the zero-temperature coupling 
of the quartic term of the bulk scalar in order to yield a convergent solution for the dilaton. 
It is noticed that as long as the parameters Cy and c\ are non-zero, even if no matter how 
small they are, the divergence at the horizon can be removed. On the other hand, as two 
parameters Cy and cx become smaller and smaller, the differences from the zero-temperature 
part given as vo{z) in Eq.l fTOl) also become smaller and smaller. Thus we will assign values 
to c^, and c\ as small as possible. 

For simplicity, we will treat Cy and c\ equally In general, they can be a 

function of temperature, but we are going to treat them as a constant also for simplicity. 



III. SPECTRAL FUNCTION FOR VECTOR AND AXIAL- VECTOR MESONS 



To begin with, we combine the gauge fields into vector field V^j and axial-vector field 



as 



V^i ^ ,:iAlM + AIm) and W^, 



2 (^i.Af 



(16) 



where M = xo,Xi,X2,X3, z as mentioned below eq.Q. The dual current operators for and 
W^.j are flavor vector current q'j^t"'q and flavor axial-vector current q'j^'j'^t°'q, respectively. 
For convenience, we will change the radial coordinate z{= 1/r) to n as 



u = z/zq = {ttT)z. 



(17) 



The equations of motion for the vector field and the axial-vector field in the 
spatial component x = Xi,X2,X2, are given as 



V : = 9, 
AV : = du 



-<S>{u) 



-$(n) 



+ 4glgUX'{u)fW^{x,u) 



+ 



^g^^d,d^V:{x,u) 



(18) 



(19) 



where "V" and "AV" denote the vector and axial- vector, respectively. As in ref.[26[, we 
adopt the axial gauge condition for the radial direction as Al,z = ^_r,2 = and the Landau 
gauge for the 4d space-time as d'^Ai,^ = O^Ar^^ = to get rid of unphysical polarization. 
Performing Fourier transformation for the four- dimensional coordinates as V^{x^u) = 

d^pe'P-^V^{p,u) and W^{x,u) = J d^p e'P-''W^{p,u) with = (oj, qi, q2, qs) and the 

momentum q taken as one direction, we have 

1 



Vr{p,u) + F,iu)V,^'{p,u) + 
W:"{p,u) + F,{u)W,fip,u)^ 



1 



Fo{u)V:{p,u) = 0, 
;Fo{u)W:ip,u)+ gl 



v?f{u) 



with 



FAu) 



a'in) ^ f'ju) 
a{u) f{u) 



$'(n), Fo{u) 



P{u) fiu) 



(20) 

W:ip,u) = 0. (21) 



(22) 



where the prime (') means a derivative with regard to u. It is seen that there is no difference 
in the solutions V^{p,u) and W^{p,u) coming from the index for SU{2) algebra a, so we 
will not attach the index a to the spectral function below. 

With the above analysis, we are now in the position to make numerical computation for 
the retarded two-point Green function G'^(a;, q) for the operators of flavor vector current or 
flavor axial- vector current. Then we can get the spectral function as 

p{uj,q) = --\mG\uj,q)d{uj^ -q^). (23) 

TT 

Here the index a will not be attached to p(a;, q) for the reason mentioned above, where 
6{(jj'^ — q^) is a step function introduced to keep the squared energy being positive in the 
boundary theory as a general requirement likewise ref. j36|. As the key point in spectral 
functions is the location and the width of the peak, where the width and the location of the 
peak are deflned in our present study by the Breit-Wigner form given in eq. fl30|) . we ignore 
the overall factor — I/tt in our actual analysis for simplicity. 

The solutions V^{q,u) and W^{q,u) described by the equations of motion fl2U|) and ([5T]) 
can be generally given by two linear independent solutions as 

V:{p,u)=A\uj,q)%{u,q,u) + B\uj,q)^l{uj,q,u). (24) 

As the axial-vector sector is the same as the vector sector, we will skip its description below. 
As described in appendix, the ratio B{(jj,q) / A{(jj,q) corresponds to the retarded two-point 
Green function, which can be written as 

^ ' A-{u,q) ^'{{u,q,u)v-{u,q,Q) ' ^ ' 

Here we have assumed that $o becomes some constant and vanishes on the boundary 
(m — )■ 0). The part "■ ■ ■ " represents irrelevant part in our present consideration. C represents 
an overall constant fixed from QCD. 

Our approach to get the Green function fl2^ is the shooting method. For that, we need 
the asymptotic behaviors for the boundary conditions. It turns out that the asymptotic 
behaviors are the same for all the models though the dilaton takes different form for the 
different models when solving the equations of motion (120|) and (!2T|) . Thus, the following 
description is common for all models. From Eqs.f l20p and (I2T!) . we can obtain the solutions 
for $Q and at the vicinity of the boundary in the limit: 2; — >■ e, with e an infinitesimal 
number. The solutions are composed of two liner-independent functions by the first-kind 
Bessel function Ji and second- kind Bessel function Fi, we assign them as follows 

V : $o(w,g,M = e) = e Fi(xv) and <l>i(cu, g, m = e) = e Ji(xv), (26) 
AV : q,u = e) = e Yi{xay) and <^1{u, q,u = e) = e Ji(xav), (27) 



_ q2 r 

with xv = — e and xav = — ^^e other hand, solving Eqs. fl20l) and 

f l2T]) around the horizon in the limit: n — )■ 1 — e, we find that the solutions for the vector 
and axial-vector mesons are the same with the so called in-going and out-going solution. As 
usual, we take the in-going solution as the boundary condition: 



V:[uj, g, n = 1 - e) = W:[uj, g, u = 1 - e) = e"^^. 



(28) 



We now make numerical computation with fixing the values for and cx and considering 
the temperature effect. In principle, one can assign their values as small as possible, while in 
the practical numerical computations, it is difficult to take extremely small values for both 
Cy = c\ and temperature T. It will be seen that the mass spectrum appears sharply like a 
spike at low temperature, while as the temperature turns on, the spike becomes wider and 
wider as well as lower and lower when the temperature gets higher and higher. In order to 
check such a behavior, we will start with the temperature as low as possible, and turn on 
temperature gradually. In the practical computations, the possible lowest temperature is 
assigned in the following figures for the different cases. In figuring out the values of c\ and 
Cy, it is seen that the positive value is excluded, which can be seen from Eq.([9]) or Eq.f lT5|l . 
for a positive value, a divergence appears in at somewhere in u moving from to 1. 

Thus, we will consider only the negative values for c\ and c„. As a consequence, the possible 
small values for ca, Cy are assigned in the practical computations to be as follows : 

The case without the quartic term (Aq = 0) : = — 10~^ (constant). ^^^^ 
The case with the quartic term (Aq = 9) : Cy = cx = — 10~^ (constant). 

Now Let us compute the spectral function with g = and compare its result with the 
data listed in Table JVl and I VII (As for g 7^ 0, see Sec|Vl) The results are shown in Figj2], 
[3] and |H As we set g = 0, the positions of the peaks in x-axis can read as the masses of 
mesons. 
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FIG. 2: The results of the spectral function for vector meson defined in eq. (|23p without the quartic 
term (Aq = 0). As g = 0, the positions of the peaks in x-axis can be read as masses of mesons. The 
parameters used here are given in eq. (|29p and Table illl and the unit of temperature is in GeV. The 
radial coordinate is taken as ti G [1 — 10~^, 10~^] in the numerical computation. The locations of 
the spike match with the results shown in Table IVl 



la axial-vector (?tQ=0) 

T=0.045 

T=0.050 
T=0.060 
T=0.075 ---- 








lb axial-vector (Xq=0) 






T=0.050 






T=0.055 - 






T=0.060 






T=0.090 









m [QeV] 



m [QeV] 



lla axial-vector (^q=0) 






T=0.045 






T=0.050 






T=0.055 






T=0.080 








''. 1 





1.0 1.5 
0) [QeV] 



lib axial-vector (Xq=0) 






T=0.046 






T=0.048 






T=0.055 






T=0.120 













1.0 
m [QeV] 



FIG. 3: The results of the spectral function for axial-vector meson. Other choices are the same 
as ones in FigjJl 
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FIG. 4: The results of the spectral function with the quartic term (Aq = 9). Other choices are 
the same as ones in Figl2]but with the parameters given Table jllll 

From the figures, it is seen that there are more spikes appearing at low temperature. In 
particular, the locations of the spikes which stand at lower temperature are very close to 
the values of the low-lying states obtained at zero temperature listed in Table IVl and IVII 
When the temperature is increased, the spikes become less and less and the peaks move to 
the low value of cu. It is noticed that the spikes collapse in turns from the higher excited 
states. Namely, the much higher a state is excited, the more unstable that state becomes. 

To check explicitly how the peaks behave when temperature becomes higher, we illustrate 
in detail the peak for the lowest lying state and plot in FigOits behavior for various temper- 
atures. It is seen that as the temperature becomes higher, the width of the peak becomes 
wider and wider and eventually the peak disappears. We shall discuss it qualitatively in 
next section. 



n pv (MeV) la lb Ila lib 



Pav (MeV) la lb Ila lib 



775.5 ± 1 739 603 777 727 

1 1465 ± 25 1223 1175 1292 1468 

2 1720 ± 20 1534 1509 1596 1744 

3 1909 ± 30 1784 1769 1842 1971 

4 2149 ± 17 2000 1990 2054 2170 

5 - 2193 2187 2249 2351 



1230 ± 40 934 714 940 807 
1647 ± 22 1468 1247 1496 1507 
1930+f[j 1822 1573 1831 1778 
2096 ± 122 2109 1829 2102 2003 



2270 



+55 
-40 



2358 2049 2338 2202 
2582 2049 2338 2202 



TABLE V: The experimental and the predicted mass spectra for vector (left side) and axial- vector 



(right side) mesons without the quartic term (Aq = 0) given in 15 1. 



n pv (MeV) Ila lib 



PAV (MeV) Ila lib 



775.5 ± 1 583 646 

1 1465 ±25 900 1468 

2 1720 ± 20 1248 1793 

3 1909 ± 30 1564 2008 

4 2149 ± 17 1860 2170 

5 2265 ± 40 2143 2289 



1230 ± 40 1128 913 
1647 ± 22 1643 1618 

1930l^[! 1953 1940 
2096 ± 122 2225 2161 

22701^^ 2486 2333 
2742 2470 



TABLE VI: The experimental and the predicted mass spectra for vector (left side) and axial-vector 
(right side) mesons with the quartic term (Aq = 9) given in jl5( |. 



From the above studies, we can draw the following statement that the peaks standing like 
a spike in the figures of the spectral function indicate the resonances of meson, which can 
explicitly be seen in Figs|2l [3]and|l]at low temperatures. When temperature increases, the 
locations of the peak move to a lower value of u and the width of the peak becomes wider, as 
explicitly seen in FigJSl Eventually, the peaks disappear completely at temperatures around 
T = 200 MeV, which is regarded as a critical temperature for the disappearance of mesons, 
or for the confinement / deconfinement transition. 

The results of FigJS] imply that the critical temperature for the restoration of chiral 
symmetry breaking in our model should also locate around T = 200 MeV. The reason is 
that the disappearance of the peak means the deconfinement of quarks within the mesons. 
In a general consideration, the restoration of chiral symmetry will occur at the temperatures 
around that point. On the other hand, the VEV v{z) is going to become much small at the 
temperatures around T = 200 MeV . This can be understood from the 5D AdS/QCD where 
the bulk coordinator z plays the role of the running energy scale, i.e., a large z corresponds to 
the low energy and a small z to the high energy. The AdS/QCD is actually more reliable to 
describe the low energy dynamics of QCD with a large z, which may explicitly be seen from 
the wave functions of the mesons as shown in Fig. 5 of Ref.[l5[, it is noticed that the peak 
of the wave function for the ground state appears at 2; ~ 1.5 GeV~^ ~ 1/667 MeV~^, when 
z < 1.0 GeV~^ the wave function goes far from the peak and becomes negligible small at 
z < 0.5 GeV~^, which corresponds to the energy scale around yU ~ 2 GeV where perturbative 
QCD gets available. Here, it is easy to check that at temperatures around T = 200 MeV, 
the allowed value for z is given hj z < l/(7rT) ~ 1/628 MeV~^ due to the requirement for 
the black hole metric / > 0. From the boundary condition at UV, it is known that with 




FIG. 5: The peak dissolves as temperature grows for the lowest lying state. The other choices are 
the same as ones in Figl3]or[31 

the small values of z, the VEV vo{z) for rriq = approaches to zero as the power z'^ for the 
given condensation a at low energies, and the part Vi{z) goes to zero exponentially. Thus 
it is expected that the VEV v{z) holographically corresponding to the chiral condensation 
with the bulk coordinator z characterizing the running energy scale approaches to be very 
small as z goes to be at the UV boundary. In this sense, we are led to the statement that 
the chiral symmetry will be restored around T = 200 MeV. 

Note that in ref. [2^ the mass spectra in vector and axial-vector have also been computed 
in a different model from ours and obtained a different dissolved temperature. We shall make 
a comparison about the dissolved temperatures between two models and discuss the physical 
reasons causing the differences. Let us first present the results for both models in Table IVlIl 
Then it can clearly be seen that two results are quite different. We now discuss the possible 
reasons causing the differences. One direct way to consider the difference is to do an analysis 
for the existence of states as in Fig.l in ref. 2^, here we will only make a general comment 



for the obvious differences in two kind of models. There are actually three differences: the 
modified background geometry due to /i^ as in eq.(j2j), the dilaton and the bulk scalar field 
as in eq.® and ([9]), and the quartic scalar field potential with the coupling A as in eq.©. 
We shall make more detail comments as follows: 

First, the introduction of the scale fi^ in the modified 5D metric will be crucial for causing 
the differences. As the differences shown in Table JVTTl are the ones at finite temperature, 
and the scale /i^ couples with temperature in the metric, which may be seen from the 
redefined new variable u = ttTz so that the modified term fi'^z'^ in 5D metric is related to 
the temperature via /i^/(7r^T^)M^. Thus tuning of /i^ term will show us how the effect of 
temperature is enhanced or softened. 

Next, turning to the dilaton and the scalar field, the dilaton in our model is determined as 
the solution of equation of motion with the given v{z) which consists of the finite temperature 



and the zero-temperature parts as in eq. (fTO!) . It is different from the case in 26| where the 
dilaton is fixed first and the VEV is determined from the equation of motion. It is important 
to notice that to cooperate the chiral symmetry breaking and linear confinement we have 
to obtain the correct boundary conditions for both the VEV at the UV boundary and the 
dilaton at the IR boundary, for this purpose we have shown in[l5| that it is necessary have 



a non-zero fig or adding a high order interaction in the scalar potential. 



Therefore, we arrive at conclusion that the deformation of IR region of bulk by the term 
^ is the key point for causing the differences given in Table JVlTl Finally, it is clear from 
Table JVTT] that the effect of the quartic scalar field potential with nonzero coupling A is the 
reason for lowering the critical temperature. 





Ao = Ao = 9 


Ao = Ao = 9 


Our results 


(lib) 


0.20 0.15 


0.20 0.14 


The results in 


ref.[26l 


0.14 


0.10 



TABLE VII: This table shows the dissolved temperatures obtained in vector (left) and axial-vector 



I— t 



sectors (right) in our present model and ref.[26||, where the unit is in GeV. Here our results are 



obtained in model lib. As the quartic scalar potential is not taken into account in ref. 26|], their 
results correspond to Aq = case. 

Before closing this section, we would like to further make comments on the several options 
considered in this paper, which include: the value for the coupling of the quartic scalar 
interaction and the choices for the regularization parameters given in eq. fl29|) . and the four 
types of models shown in tableHl First, as for the coupling of the quartic scalar potential, 
as mentioned in the introduction, it was first introduced for obtaining a nontrivial solution 
of dilaton and meanwhile yielding the chiral symmetry breaking, while it was pointed out 



that the bulk scalar potential field is not bounded due to the opposite sigen |14l] . Concerning 



the values assigned in this study, we have followed the analysis in ref. [15j with changing the 
sign for the stability consideration and the nontrivial solution of dilaton is obtained as the 
consequence of the modified 5D metric at IR region with the introduction of the yU.^ term, 
and the value of the coupling A is taken from ref. 1^ by fitting to the experimental data. 
Concerning the regularization parameters Cy and c\ as shown above, we have assigned the 
values as small as possible in the numerical calculations, so that they do not affect our final 
conclusions. As for the four types of models discussed originally in ref. [l5| , which is to show 
how the predicted mass spectra rely on the IR boundary behaviors and the different forms 



with the same boundary behavior of the VEV. As a consequence, it can be seen from ref. [15 
that the results are more sensitive to the boundary behaviors classified as 'a' and 'b' than 
the different forms named as T and 'IF in TableJH the model lib leads to more consistent 
results with experiment data than other types of models. This is the reason why we take 
the model lib as the studying model in this paper. 



IV. MASS SHIFT AND WIDTH OF PEAK FROM SPECTRAL FUNCTION 



As we show in previous section how the peak dissolves as the temperature is increased 
to a critical point. Here we are going to make a numerical study with fitting the spectral 



function by the following Breit-Wigner form: 



au 



(a;2 - m2)2 + 



(30) 



Here m and F correspond to the location and width of the peak. P{uj'^) represents the part 
not forming peaks and is taken to be the form Pioj"^) = Ci + C2U}'^ + c^iuj^Y'^ ■ Thus it involves 
eight fitting parameters: a, 6, m, F, Ci, C2, C3 and C4. To see quantitatively the mass shift 
and the width change of the peak, we consider the lowest lying state and plot the relation 
between the mass shift Am^ = (m — mo)^ and the width F against the temperature. Here 
mo is the predicted meson mass at zero temperature given in Table fVl and I VII The result is 
shown in Figsj6] and [71 It is seen from Figj6] that Am^ is proportional to F as shown by the 
dotted line. 

From FiglSl it can explicitly be seen that the mass shift Am^ grows linearly as temperature 
increases and begins to deviate from the linearity near the temperature which leads the 
spectral function to dissolve. A similar phenomenology was also observed in ref. 26| and 



considered to be consistent with the results based on QCD sum rule calculation [371 . 
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FIG. 6: The relation between the mass shift Am? = (m — mo)^ and the width F, which is obtained 
from the fitting to the numerical results of spectral function for the lowest lying state by the Breit- 
Wigner form given in Eq. ([30|) . The parameters used here are given in Eq.(|29p and Table HIl and 
mH and the unit of temperature is in GeV. The dotted lines with the linear behavior are plotted 
for a comparison. 
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FIG. 7: It shows the temperature effect in Am = m — mo and F. The other choices are the same 
as the ones in Fig|6l 

It is noticed that the relation between Am^ = (m — mo)^ and F is not known in the low 
temperature region where a peak appears as a spike. This is because there is no numerical 
results as in such a low temperature region, the peak stands like delta function, it is difficult 



to measure the width numerically. From FiglTJ it is seen that the linearity occurs above 
certain temperature. In the low temperature region, such a linearity changes since the mass 
m and width F get values m = rriQ and F = at zero temperature, respectively. Thus, it is 
expected that a similar change happens in the relation between Am^ = (m — mo)^ and F. 



V. FINITE MOMENTUM EFFECT 



In this section, we are going to study the finite momentum effect in the spectral function. 
Such a momentum effect in the shape of spectral function was also studied in ref. [26|, 27 . 
Our results are shown in FigJHl Note that the spectral function in our present consideration 
is defined with a step function in Eq.f l23p for a positive squared energy in the boundary 
theory from a general point of view (36^ 
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FIG. 8: The momentum effect in the spectral function for various momentum in unit GeV with 
T = 0.10 GeV. The parameters in the left and the right are given in Eq. (j29p . Tablellll and lllli The 
radial coordinate are taken to be n G [1 — 10~^, 10~^] in the computation. 

It is seen that there is a wall-like behavior caused by the step function defined in Eq.( l23|l . 
In particular, when the momentum turns on, the peak moves to large direction and 
gradually decays. Such a phenomenology is consistent with the observation in jsM 



39|, but 

different from the results obtained by the perturvative analysis in a field theory [40|, which 
may be considered as a peculiar and common behavior in strongly coupled field theories. 

Let us turn to the effect of the step function defined in Eq. fl2^ . Its effect appears as a 
wall in FigJHl It is seen that there is a blow-up around the wall, and the hight of the wall 
is changed depending on fineness of step of taken in the numerical computation. As a 
result for the introduction of the step function, one may notice the difference in the speed 
of the wall moving as is increased. It may be seen from the differences between the left 
plot and the right plot in Figj8]with Aq = and Aq = 9, as well as with parameters given in 
TableJUandlml 



VI. CONCLUSIONS AND REMARKS 



In this paper, we have extended the predictive soft- wall AdS/QCD model [15|, |l6| to a 
thermodynamic model by introducing a black hole metric. As the predictive AdS/QCD 
model has adopted an IR-improved 5D metric to characterize the low energy dynamics of 



QCD [15|, |16|, we have shown that it has to figure out a modified bulk vacuum expectation 



value and a modified coefficient of quartic term for the scalar field to get a non-singular 
dilaton solution. In obtaining the modified bulk vacuum expectation value and a modified 
coefficient of quartic term, an additional term is introduced to regulate the singular behavior. 
For that, the coefficients of the additional terms are taken to be small as much as possible in 
the practical calculations. With the smooth dilaton solution, we have computed the spectral 
function for the vector and axial-vector mesons. 

It has been demonstrated that the peaks in the spectral function indicate the resonance 
mesons, as shown in FigsEj [3] and HI The location and the width of the peaks moves to a 
lower value of u and becomes wider as the temperature increases, respectively, as seen in 
FiglSl We have found that the peak disappears completely at temperatures around T = 200 
MeV, which implies the deconfinement of quark and the restoration of chiral symmetry 
breaking. To see how the peak dissolves quantitatively when the temperature is increasing 
to the critical point, we have made a numerical study by ffiting the spectral function with 
taking the Breit-Wigner form. It has been seen that the mass shift Am^ = (m — mo)^ with 
rriQ as the meson mass at zero temperature grows in a linear way as temperature increases. 
We have also investigated momentum effects and shown that the peak moves to a large u"^ 
direction. 
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Appendix A: The formula for two-point retarded Green function 

In this appendix, we will show how to obtain Eq.f l25p . It is seen that the Lagrangian up 
to third-order perturbations of vector and axial-vector sector can be written as 

W = - Jd'pe'^-S^f{z){K\p,z)Y + j^^{t^{p,z)Y^ + ..., (Al) 

where "■ ■ ■ " denotes the axial-vector part which has the same form as the vector part and 
we abbreviate it. 

In general, the action is assumed to have the following form 

S[<P] = jdz {^(^)0(^)0"(^) + i3(^)(0'(^))2 + C(^)0(z)0'(^)+P(^)(0(^))2}, (A2) 

where (t){z), A{z), B{z), C{z) and T>{z) are arbitrary coefficients. It can be shown that 
Eq. (lA2p may be rewritten as follows 

= 1 [dz |eOM ■ (f){z) + 9,(ST)}. (A3) 



"EOM" denotes equations of motion and it will vanish as the bulk gravity is evaluated at 
classical level. On the other hand, "ST" denotes a surface term given as 

ST = A{z)(f){z)(l)'{z) + {(l){z)A{z) - (l){z)A'iz)}(l)iz) + 2S(z)0(z)0'(2) + C{z)(l)\z). (A4) 

Here we will mention briefly how to obtain eg. (1A3I) from eg. (1A2I) : Firstly, we take the 
variation in an usual way to get eguations of motion as S[(j) + S(f)] — S[(j)]. Then it can be 
summed up in "EOM" part and "ST" part up to the second order, i.e., S[(j) + 6(f)] — S[(f)] = 

J dz I EOM ■ 50(2) + 9^(ST)| + 0{S(f)'^). Then, discarding the higher-order terms and 

replacing 6(f) with (f>, its Ihs becomes 2S[(f>] (provided that powers of all fields are same), and 
we arrive at Eg. (lA3p . 

Using this prescription, eg.( 1Aip can be written at classical level as 



-''^^^a{z)f{z) 



augc 



d^pe'P-''V:{p,z)Vf{p, z) 



The evaluation at ^ = in the above can be written as 

e-*(^)a(;2)/(2) 



+ ■ 



z=0 



d^p e'P-''V^{p,z)V,^'{p,z) 



2=0 



,4 ^r>■. V.\P,z)V:'iP,^) 

d p 



z=0 



d^p e'P 



,^A%u,q)A'^{-u,-q) 



9l 



_A-{u,q) 



2 a;2 



(A5) 



7£; + log -Va;2 - g2 



(A6) 



z=0 



Here we have substituted eg.( l26i) . A'^{uj,q) = A"-{—uj, —q) as can be seen from eg.( fT8|) and 
( 12T|) . '-jE means Euler constant. Thus following GKP-W relation and the prescription given 

in ref.jij, the retarded Green function = j e-''-^^(t)( [J^(x), J^(0)] ) (As 

for what J°(x) means, see below eg.l lTB]) .) can be obtained as 



B'^{u,q) ^q 



2 a;2 



gl [A-{u,q) 



|7s + log _ I 



(A7) 



z=0 



Thus from the imaginary part given in B"-{uj,q)/A°'{u,q), we obtain the spectral function 
for vector and axial-vector mesons. 
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